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1 Introduction 



Given some continuous differential equation, we want to find its discretizations 
which preserve some global property, namely the explicit linearizability or more 
generally the discrete Painleve property, a notion recalled section 2| together with 
its group of invariance. From this point of view, we will never have to address 
the question of finding a continuum limit, and consequently we will always write 
the discrete equations in a form as close as possible to the canonical form of the 
continuum limit, which is well established [^, 20[| . 

Discrete equations are functional equations linking the values taken by some field 
variable m at a finite number iV + 1 of points, either arithmetically consecutive : 
X + kh, or geometrically consecutive : xq^ , k — = 0,1, . . . , N, where or g is the 
lattice stepsize, assumed to lie in some neighborhood of, respectively, or 1, and k^ 
is just some convenient origin. The integer is called the order of the equation, 
and we denote, for brevity, d-(E) and q-(E) these "difference" and "g— difference" 
equations. Their study was initiated by Laguerre, mainly as three-term [N = 2) 
recurrence relations between coefficients of orthogonal polynomials. This remained 



for long a mathematical subject ||5^, |T9|, which then extended to topological field 



theory 0, . Finally, the discrete equation 

E = -{u-2u + u)/h'^ + 2{u + u + u)u + x = (1) 

already considered by the authors of last five references, was again encountered by 
statistical physicists in two-dimensional quantum gravity P, who recognized 

it as a discrete analogue of the first Painleve equation (Plj 

E = -u" + + X = 0. (2) 

The same happened simultaneously with a discrete analogue of the second Painleve 
equation (P2) Q (in the particular case a = 0) 



E = -{u-2u + u}/h^ + {u + u)u^ + xu + a = (3) 
E = -u" + 2u^ + xu + a = 0. (4) 

The above short notation 

even : u = u{x), u = u{x + h), u = u{x — h), u = u{x + 2h), ... (5) 
odd : u = u{x + h/2), u = u{x - h/2), W = u{x + 3h/2), . . . (6) 

is adopted throughout the article. 

The reason why these two discrete equations, among many others with the same 
continuum limit, deserve the name of discrete Painleve equations, in short d-(Pl) 
and d-(P2), is that they possess the discrete Painleve property. Indeed, both admit 
a discrete Lax pair. 

Up to now, there exist two methods to find discrete Lax pairs : the discrete 
isomonodromic deformation method 0, and the discrete analogue of the method 
of Zakharov-Shabat and Ablowitz-Kaup-Newell-Segur (AKNS). But both methods 
are inverse methods, i. e. they generate some discrete equation as a condition betwen 
the coefficients of two given linear operators. The drawback is that the obtained 
discrete equation may not be the one which was looked for. We propose here a direct 
method, based on discretization rules, to search for the Lax pair of a given discrete 
equation, and we obtain several new Lax pairs in this way. 

Just like its continuous counterpart, the discrete Painleve test is the set of all 
methods one can imagine to build necessary conditions for a given discrete equation 
to possess the discrete Painleve property. Two such methods are known, the sin- 



gularity confinement method ||2J] and the method of perturbation of the continuum 
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limit [O. The latter is based on the perturbation theorem of Poincare, which is 
applicable to differential systems in an arbitrary number of independent variables, 
whether discrete, continuous or even mixed discrete-continuous. 

The paper is organized as follows. The discrete Painleve property (PP) and its 
group of invariance are defined in section ^, and a first set of basic rules of discretiza- 
tion is given in section ^. Section ^ recalls as an illustration the exact discretizations 
of the elliptic equations. In section ^, we define discrete Lax pairs from continuous 
Lax pairs. Section |^ states necessary discretization rules for Lax pairs and details 
the direct method to obtain the discrete Lax pair of a given equation, with an appli- 
cation to the d-(Pl) (0), the d-(P2) (Rj) and a particular d-(P3) of degree one. In 
section ^, this method is extended to ciiscrete equations with complementary terms 
which do not contribute to the continuum limit, and several new Lax pairs are ob- 
tained. Section R] explains the method of perturbation of the continuum limit for 
the discrete Painleve test, and applies it to a qualitative candidate d-(Pl). Finally, 
we define in section |^ criteria for a systematic search for the fifty discrete Gambler 
equations. 



2 The discrete Painleve property and its group of 
invariance 

The (continuous) Painleve property is defined as the absence of movable 
critical singularities in the general solution of a differential equation 

Vx : E{x,u,u',...,u^^^) = 0, (7) 

where a singularity is said movable (as opposed to fixed) if its location in the complex 
plane of x depends on the initial conditions, and critical if the solution is multivalued 
around it. For shortness, following Bureau 0|, we will use the terms "stability" for 
PP, "stable" or "unstable" for an equation with or without the PP. 
The PP is invariant under the group of birational transformations 

{u, x) (f/, X) : u = r(x, U, dU/dX, . . . , d^'^U/dX^-^) = 0, x = H(X),(8) 
{U, X) ^{u,x): U = R{X, u, du/dx, . . . , d^-\/dx^-^) = 0, X = ^(x), (9) 

(r and R rational in U,u and their derivatives, analytic in x,X). An easier to 
manage subgroup is made of the homographic transformations 

(u, x) ^(U,X): u = ^J^^, X = ^(x), ad-bc^O (10) 

cU + d 

where (a, b, c,d,C,) are arbitrary analytic functions of x. In his classification of second 
order first degree equations. Gambler ||20| has found respectively twenty-four and 



fifty equivalence classes for these two groups (with minor later corrections p, ) 
In the discrete case, let us consider equations 

Vx Wh : E{x, h, {u{x + kh), k - h = 0, . . . , N}) = (11) 
Vx Vg : E{x, q, {u{xq''), k - ko = 0, . . . , N}) = (12) 

algebraic in the values of the field variable, with coefficients analytic in x and the 
stepsize h or q. It should be noted that m is a function of two variables, x and the 
stepsize. A natural definition for the discrete Painleve property is the following |l^ . 



Definition. A discrete equation is said to possess the discrete Painleve property 
if and only if there exists a neighborhood of h = (resp. g = 1) at every point of 
which the general solution x u{x,h) (resp. x u{x,q)) has no movable critical 
singularities. 
Remarks. 



4 



1. The definition reduces to that of the continuous PP in the continuum limit. 



2. The singularities in the definition belong to the complex plane of x, not of the 
stepsize. 

3. This definition immediately extends to equations in an arbitrary number of 
independent variables, discrete or continuous, the extension starting then from 
the definition of the PP suited to partial differential equations (PDEs), which 
we do not remind here since this is not our subject. 

The discrete PP is invariant under the discrete analogue of (P), which is the 
group of nonlocal discrete birational transformations 

u = r{x, h or q, U, U,U_, . . .), 

U = R{X, H or Q,u,u,u, .. .), X = C,{x, h or q), H = rj{h), Q = K{q), (13) 

(r and R rational in U,U,U_, . . . ,u,u,u, . . ., analytic in x and the stepsize, ^, rj, k 
analytic). There exist two discrete analogues of the subgroup (|l3), and both may 
be useful to establish the discrete equivalent of the classification of Gambler. The 
first one is the group of transformations (|T^) which in the continuum limit reduce 
to the homographic transformations (0), where {a,b,c,d,^) are arbitrary analytic 
functions of x and of the stepsize. The second one is the group of local homographic 
transformations (r and R homographic in U and u, independent ofU,lJ_, . . . ,u,u, . . ., 
analytic in x and the stepsize, ^, r], k analytic). 
Remarks. 

1. The first subgroup seems more useful, although it does not contain the trans- 
formation u = h''U, k & Z. 

2. Just like in the continuous case, the birationality can only be proven by taking 
the discrete equation into account. For instance []TB[, given the equation 

{u + u){u + u){4u + h^x-6) + 8 = 0, (14) 

and the transformation 

U{X) = 2/{u{x-h/2)+u{x + h/2)), (15) 

one first deduces the inverse transformation 

u{x) = (6 - 2U{X - H/2)U{X + H/2) - H^X)/A, X = x, H = ^16) 

then one plugs it into the direct transformation to get the transformed equation 

{U + U)U'^ + {H*X -6)U + 4 = 0. (17) 

The fields which admit a continuum limit are {u — l)/h^ and (1 — U)/H'^, this 
limit being (PI) for both fields. 



3 Basic rules of discretization 

Let us now give some basic rules for discretizing a given continuous equation (|^) 
into either a difference equation (|ll|) or a g-difference equation (|12D. 

The question of discretization is well known in numerical analysis, where one 
looks for a scheme of discretization which maximizes the order, called scheme order, 
of the remainder of the expansion of the left-hand side of (0) in a Taylor series of h 
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around the center of the N + 1 points. A scheme of discretization is said exact iff it 
has an infinite order, hke that for the particular (P3) equation a = (3 = 'y = S = 

— uu" + u'"^ — uu /x = 0, u = CiX^"^ , u{xq)u{x / q) — u{xY = 0. (18) 

Contrary to numerical analysis, only interested in a local integration, we require 
the scheme of discretization to preserve the differential order N, an essential element 
for a global knowledge of the solution : every discretization must involve + 1 
consecutive points. 

Like in the continuous case, there exists another important element concerning 
discretization rules. 

Definition. The degree of a discrete equation is the highest of the two polynomial 
degrees of the LHS E of the equation in u{x) and u{x + Nh), or u{x) and u{xq^), 
where E is assumed polynomial in the iV + 1 variables u{. . .). 

The following conjecture has recently been made |jl3[ : "Given an algebraic 
differential equation with the PP, there exists a discretization scheme of order two 
which conserves the degree. " 

This conjecture was supported by two examples, a second order first degree 
equation with the PP 

u = (cix + C2)^ uu" - (1/2)m'2 = 0, (19) 
and a first order first degree equation without the PP 

u = {x- Cl)-^/^ u' + {l/2)u^ = 0. (20) 

Both admit an exact algebraic discretization, resulting from the elimination of Ci 
or (ci, C2) between the values of u taken at two or three contiguous points, and the 
resulting discrete equation has degree two. For the first example, we constructed a 
first degree discrete equation with the PP, while we showed the impossibility to do 
that for the second equation. Another example is displayed in section 

This led us to state the additional rule, restricted to algebraic differential equa- 
tions with the PP : every second order scheme must also conserve the degree. 



4 The exact discrete elliptic equation 

The (continuous) elliptic equation has two usual kinds of normalized forms, the 
one of Weierstrass and the twelve ones of Jacobi, defined by the first and second 
order equations {g2,g3,a,b constants, (p,q,r) arbitrary permutation of (c,d,n)) 

p'2 = 4p5-g^p-g^ (21) 

p" = Qp'-g2/2 (22) 
ps'^ = qs^rs^ = (ps2+a)(ps^+6) (23) 
ps" = ps(qs^ + rs^) = ps(2 ps^ +a + b). (24) 

Elliptic functions possess an addition formula, i. e. an algebraic relation between the 
values of the function at the three points {xi,X2,xi + X2), with (xi, 3:2) arbitrary. 
As noticed in the context of discretization by Baxter [§] and Potts [Q, the choice 
{xi, X2) = {x — h/2, h) ipso facto defines an exact discretization scheme for the first 



order equations (|2l|) and (pSj) . A scheme for the second order equations ( ^21) and 



41) then results nx)m the difference of the discrete first order equations taken for 
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h) and {x — h,h). The results are |4^, with notation 
751) and (|23), (|) for (|6|) and (@) 



(u — 2u + u)p{h) 
{u — 2u + u) ps^{h) 



2uu{u + u) - (5-2/2) (m + u) - gs 
-[{uu + 02/ 4f + 93(11 +u)]p-\h) 
2u{u + u + u) — (72/2 
-[u\u + u) + {g2/2)u + g3]p'\h) 
(uuy - 2(ps'(/i) + ps^{h))uu + ah 
u^(u + u) - 2(ps'(/i) + ps^{h))u 



for 



(25) 

(26) 
(27) 
(28) 



Remarks. 



1. The general solution of (|25|) and (|27|) is by construction p{x — Xq, 92,93) and 
ps(x — xo,k), where the step h is arbitrary, i. e. not necessarily small, with k 
the Jacobi modulus. These equations therefore possess the discrete PP. The 
equations (^) and (p8|) also possess the discrete PP since they admit a discrete 
Lax pair, see (|6^)-([68[) and (|62D-(|6^). 



2. Order and degree are conserved by the four discretizations. 

3. From these schemes of infinite order, the Laurent expansion of p{h) around 
its double pole /i = 0, or of ps(/i) around its simple pole, defines the second 
order schemes 



[u - 2u + u)h ^ 
{u-2u + u}h'^ 

which have the discrete PP, see 



2u{u + u + u) — (72/2 
u^ivi + u) — {a + h)u, 

-(0) and (§2|)-(^. 



(29) 
(30) 



4. Similar equations hold for the nine other Jacobi functions pq, where the coef- 
ficients of the r.h.s. of ( pTj ) and (p8|) are polynomials in /c^, pq^(/i), pq'(/i). 

5. Equation (P^) was obtained and integrated in 1973 by Baxter in the eight- 
vertex model, as a commutation condition of the Yang-Baxter, or star-triangle 
relations. These Yang-Baxter relations PS[, which are second order discrete 
tensorial equations, play in the discrete oomain a role as central as the one 
played by the Yang-Mills equations in the continuous domain. 



5 Discrete Lax pairs 

In the two relations defining a Lax pair of a (continuous) ODE E{x,u) = 0, 
e.g. in matrix form with t the spectral parameter 

d^ip = Li), dtip = Mip, (C = dtL - d^M + LM - ML = 0) ^ (E = 0)(31) 

one can discretize either x alone or x and t. Let us restrict here to the first case and 
to the difference type equations. 

The rule of conservation of the differential order requires the column vector 
to be discretized with two points, which we denote ip = ip{x + h/2) and ip = 

il){x — h/2), and L to be discretized with as many points u as required by the 
differential order of the equation under consideration, points which we denote u = 
u{x + h),u = u{x),u = u{x — h) for a second order equation. In order to keep a 
linear correspondence between the continuous operators (L, M) and their discrete 
counterparts, it is then convenient to discretize dxip = Lip in a dissymmetric-looking 
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way and to introduce the linear operator A linking ip to ip, thus defining the 
discrete Lax pair {A, B, z, ip, h) as ||26 



= A'4)_, = Bf_, {K = d,A + AB - BA = 0) ^ {E = 0). (32) 

The continuum limit is then 
A — 1 

— >L, (dz/dt)B^M, 

h 

{dz/dt){d^A + AB- BA)/h dtL - d^M + LM - ML, (33) 

with some link F{t, z,h) = between the spectral parameters t and z. 

For a second order equation E{u,u,u,x,h) = 0, the operators A and B must 
have the m— dependences A{u,u,u), B{u,u). 

Remark. The definition is invariant under the involution 

{E, A, B, X, h,u, u, u) {E, A~'^, B, x, —h,u, u, u), (34) 

which sometimes allows to suppress an undesired denominator in matrix A, like in 
the matricial Lax pair of d-(Pl) given in Ref. [pT] . 

The interest of a discrete Lax pair is to provide a constructive proof of the PP, 
just like in the continuous case. 



6 A direct method towards matricial Lax pairs 

The two methods recalled in the introduction to find discrete Lax pairs are 
inverse methods. From the point of view of discretization, a new, direct method 
emerges, which is as follows. 

Let be given a continuous equation, its Lax pair, and some discretization of 
the continuous equation. We first state general rules for discretizing the Lax pair, 
involving some free functions; we then enforce the cross-derivative condition (|32) to 
remove all the freedom on these functions, except the link between t and z which 
cannot be removed; finally, we choose the link between t and z, so as to preserve 
the rational dependence of the matricial Lax pair on the spectral parameter and to 
ensure the existence of the continuum limit of the discrete Lax pair. 

In this section, we handle the d-(Pl) and the d-(P2) (|^), whose matricial Lax 
pairs with a continuum limit have been found respectively by Refs. using 



an inverse method, and by Ref. [|l^], using the present method. We also handle the 



d-(P3) of Ref. [^] when its degree two reduces to one, i.e. for (a, (3, 7, 6) = (0, 0, 0, 0) 

E = -xuiu -2u + u) /h^ + x{u -u){u- u) /h^ - u{u - u) /{2h) = 0, (35) 

with continuum limit either of the two canonical forms of the third Painleve equation, 
(P3) or (P3') (Ref. m p. 1115) 



(P3')i^ - -u" + ^-^ + ^^^l-^ + f + l = 0, (36) 
u X Ax Au 

,^ , ^ ,, u'"^ u' av? + B n 5 , _, 

P3) E = ~u" + + + 7^3 + _ = 0. 37 

U X X u 

One must first make a choice between three kinds of Lax pairs for the (Pn) 



equations : the second order scalar "Lax" pairs of Garnier the second order 



matricial ones of Jimbo and Miwa ||2^, and the ones of Flaschka and Newell [jl 
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arising from the reduction of a PDE. In this paper, we restrict to the third type in 
matricial form when the matrix order is two. Using Pauh matrices 

^1 = (J J) , ^2 = (^^ , (T3 = Q , a^ak = 5jk + iejki(Ti, (38) 

1\ _ /O 



^ =' Oj' ^ = VI 



the continuous pairs are as follows. 
For (PI) : 



L = 2{u-t)a+ + a-, (39) 
M = 2u'(r3 + (-4u^ -2x + 8tu - 16t^)(r+ + (4n + 8t)a", C = 2a^E, (40) 



For (P2) |T 



L = -t(T3 +MCri, (41) 

M = (4)^2 - 2u^ - x)a-i + 2u'ia2 - (4tM + at~^)ai, C = 2ia2E, (42) 
For the (P3) : 

L = {l/2){u'/u + cu + d/u)(T3 + tai, = 'y,cP = -S, (43) 
M = {x/{2t)){u'/u + cu + d/u)a3 + xcri+r^cdxai/2-t~^{ac- Pd)a3 
—t^'^{au + '-yxu^ + c{xu' + M))cr~*'/2 

-t-'^ilSlu + + c/(x(1/m)' + l/u))a-/2 (44) 

+ cd 

C = {{x/{2tu)){a3-{cu/t)a+ + {d/{tu))a-)-^^^E. (45) 

This Lax pair for (P3) is the extension to arbitrary a, (3, 7, 5 of the pair given by 
Milne g for 7^ = 0. 

The rules of discretization are the following. 

1. conserve the matricial order. This is indeed the differential order of the scalar 
Lax pair, which must be conserved; 

2. replace the continuous spectral parameter t by an unspecified function T{z, h); 

3. discretize the operator L centered at the three points x — h,x,x + h. If L is 
traceless, so is its discretization; 

4. discretize the operator M centered at the two points x — h,x. If M is traceless, 
so is its discretization; 

5. for each matrix element, enforce conservation of order and degree; 

6. replace each monomial (du/dx)'^ by its discretization obeying the general rules, 
multiplied by the k-th power of an unspecified function g{z, h). This function 
g, whose continuum limit must be 1 for any z, represents the ratio of the 
stepsize h to the differential element dx; 

7. take B as the product of the discretized M by an unspecified function J{z, h) 
(like Jacobian) representing a discretization of the derivative dT/dz; 
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8. take {A — l)/h equal to the sum of the discretized operator L and a diago- 
nal matrix of unspecified functions of {z, h) only, dia.g{gi, g2); these functions, 
whose continuum limit must be zero, account for the dissymmetry of the for- 
mula defining A. 

In the second order, first degree case of Painleve equations (Pn), examples of 
discretizations obeying the above rules are 

in L Piuiu + m)/2 + + z/smm, i/i + z/2 + z/3 = 1, 
in M — > uu 

X in M X- h/2. (46) 

Finally, after the cross- derivative condition has been enforced, one must perform 
the continuum limit on those of the functions T, g, J, (71, (72 which remain unspecified. 
One such free function is T, because the choice of z is arbitrary. One chooses T(z, h) 
as a rational function of z (to conserve the rational dependence of the Lax pair on 
its spectral parameter) such that the inverse function z of (t, h) admits for every t 
a finite nonzero limit when h ^ Q. 

The discretized matricial Lax pair of the d-(Pl) (|l]) is, with Ai + A2 = 1 



A = 1 + h 



gi 2A1M + A2(m + m) -2T 

1 92 



R/7 _ ( 2g{u-u)/h -Auu-2{x-h/2)+4:T{u + u)-16T^ . , . 
' ~ y2{u + u) + 8T -2g{u-u)/h ^' 



and that of the d-(P2) (D is 

A = 1 + h{-Ta3 + {\iu + \2iu + u)/2)ai) + h 



91 
92 



B/J = {-2uu + AT'^ -{x-h/2))a3 + 2g{{u-u)/h){i(X2) 

-{2T{u + u) + a/T)ai (48) 

depending on the yet unspecified functions T, g, J, gi, g2 of (z, h). The cross-derivative 
condition is enforced by eliminating for instance u (or here the variable a;, which 
appears always at the first power) between the condition K = Q (^) and the dis- 
crete equation, and identifying to zero the resulting polynomial in the three variables 
(x, u,tt) (or here the three variables (m, u,tt)). This results in the relations 

d-(Pl) Ai = 1^ g, = g2 = {g- 1-9^- 2h'T' = 0, J = -g'/h', (49) 

d-(P2) Ai = 1, g, = g2 = {g- l)/h, 1 - / + h'T^ = 0, J = g'/i^T), (50) 

and one of the five functions remains free, as expected. A convenient choice of T is 

d-(Pl) T = z~ h^zy2, z = t + {t'/2)h^ + 0{h^), (51) 

d-(P2) T= , z = l + ht + Oih^). (52) 

2h 

Finally, the matricial Lax pair of the d-(Pl) (|1]) is 
l-h'^z h{2u-2t)\ 22 



A = 'Y_,^2, ), t = z-h'z'/2, z = t+{t'/2)h' + 0{h') 

f2{l-h'^z){u-u)/h -4uu-2{x - h/2) +At{u + u) -lQt^\ 

[ 2{u + u) + 8t ~2{l -h^z){u-u)/h j'^^^^ 

h~^K = 2a-iE, 
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a result obtained [O] with this method. This is the rewriting of eqns. (1.12)-(1.13) 
of Ref. [p3] which clearly shows the continuum limit, under the change of basis 



(54) 



C = (1 - h^z)/h, uj = {l- 2h\)/h^, l/f = Vl - 2h^u. 
As to the matricial Lax pair of the d-(P2) (^) 

hzB = {-2uu + 4t'^ -{x-h/2))cx3 + {z + l/z){{u-u)/h){ia2) 

-{2t{u + u) + a/t)ai (55) 
zK = 2ia2E. 

it was already known p2| . 

The case of d-(P3) is slightly more difficult due to denominators u and x in 
the expression of u", so the discretization of u'/u in L may be anything in between 
discr(M')/ discr(ti) and discr(a;uu')/ discr(a;ti^), where discr() symbolizes a discretiza- 
tion obeying the rules, and the same also applies to the term xu'/ u in M. This dif- 
ficulty is overcome by noticing that the only term of the commutator C to generate 
the equation is M^, due to the first integral xu' /u of this very particular (P3). This 
first integral has the discrete analogue /(x) = {2/h)x{u{x+h/2)—u{x — h/2))/{u{x+ 
h/2) +u{x — h/2)), so the discretization of xu /u in M can only be g{z)f{x — h/2) 
according to our rules. We therefore start from the assumption 

A = (^n n ]+h{gF{x)a,/2 + Ta,), 



G2. 

B/J = g{f{x-h/2) + fo)a3/{2T) + {x-h/2)a,, (56) 

with six unknown functions of z (T, g, J, Gi, G2, fo) and one unknown function of x 
{F) representing the discretization of u'/u. The condition that the discrete equation 
be a factor of the commutator K (|3^) yields 

d-(P3) Gi = G2, = 0, = 0, J = G'Jih'T), {Gl - h'T^y = 0, 

F{x) = {f ix + h/2) + f{x - h/2) + 2Kq)/{2x), Kq constant, (57) 

resulting in the same convenient choice of T{z) as in d-(P2), eq. (^2]). 
The resulting Lax pair for this particular d-(P3) is finally 

. ^ u{x + h/2) - u{x - h/2) 2 ^ f(^x + h/2) + f{x-h/2) + 2Ko 

•^^^^ ''u{x + h/2)+u{x-h/2)h'^''' 2x 
A = {z + z'')/2 + h{Fix)a,/2 + iiz-z-')/i2h))a,), 
h 

hzB = - — —f{x-h/2)a3 + {x-h/2)ai, 

zK = -^—^{{z + z-^)as/2 + hf{x)ai)E. (58) 
Remarks. 
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1. The scalar Lax pair satisfied by the second component il)2 of in ( p3D 

- 2^)2 + ±2)/h^ - 2Mth + 2^V^2 = 0, 

1p2,z + 2[Az{i>2 -i>2 + 2uil)2 - U^p2 - mp2]/h 

+ 4[z\^2 - 3^2) - zutlj2]h + Az'^^2h^ = (59) 

is of course the one obtained by Joshi et al. [M; as to the scalar Lax pair 
found in quantum gravity |2^, which contains algebraic coefficients, it 
is the transformed of this one by ■i/'s = Gip2, with G satisfying the discrete 
equation G/G = \/l — 2h'^u. 

2. The term u' / u of operator L for (P3) is discretized as the quotient of discr (xmm') 
by discr(xM^), a result not easy to guess in advance. 

3. The number of singular points in the complex plane of the spectral parameter 
is not necessarily the same for the matrix M{t) and for its discretized B{z). 
For (P2), M has two singular points t = 0, cxo while B has four such points 
z = 0, 00, 1, —1. 

One similarly obtains a Lax pair for the discrete Weierstrass equation without 
complementary terms pPj ) 

i = (l-AV(2/i=), 

•4 = (ft '"r")- 

B = 2X{{u-u)/h)a3 + {-m^ + At{u + u) + g2- 4:uu)a^ 

+ {2{u + u) + 8t)a' (61) 
h-^K = 2a^E 

and one for the discrete Jacobi equation (pO|) 

t = (A-l/A)/(2/i), 

' - it T). («^) 

B = {-2uu + At'^-a-b)(T3 + {X + l/X){{u-u)/h){ia2) 

-2t{u + u)cri (63) 

h-^K = 2ia2E. 

These last two Lax pairs do not depend on the spectral parameter z, and A is 
an arbitrary constant. 



7 More on discrete Lax pairs 

Some discrete equations have complementary terms which do not contribute to 
the continuum limit, such as the discrete Weierstrass equation with its terms p~'^{h\ 
or the following d-(Pl) (Ref. |18| eq. (2.8)), which only differs from the d-(Pl) (|l|) 
by terms homogeneous to K^xu and h'^u^ 

E = -{u-2u + u)/h'^ + 2{u + u + u)u + x 

-h^xu- h^u^(u + u + u) = 0. (64) 
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A Lax pair has already been obtained |]T8|, not for this d-(Pl) exactly, but only 
for the discrete derivative of a birational transform of it, which makes it not simple 
at all. Let us obtain the natural Lax pair. 

One assumes the same qualitative form ( ^Tj) than for the d-(Pl) without com- 
plementary terms, and one adds to A and B as many matrices of order at least h as 
there exist divisors of these complementary terms, each matrix being the product 
of such a divisor by a matrix of free functions of z. The result for is 



t = {1- z)/h\ z = l- h% 

^ ~ [{h/2){z + l-h^u) z j' ^^^^ 

-h^zB = 2z{{u-u)/h)cy-i-h{2uu + x-h/2)(^^^^ ^/^^ 

h-^K = 2a^E. 

This Lax pair admits by construction as its continuum limit. 

One similarly obtains a Lax pair for the discrete Weierstrass equation with com- 
plementary terms (^) 

t = {l-\)/H\ H-' = p{h), 

A _ [ A {2/H){X-l+H'u)\ 

^ ~ \{H/2){X + 1- H^u) X j' 

B = 2X{{u-u)/H)as-H{2uu-g2/2-H'gs)(^^^^ 

-2t((l + A)/A - {H% + H'g,)lA) ( ^ 
H-^K = 2a^E. 



8 The discrete Painleve test 

In the continuous case, all the methods of the Painleve test, without exception, 
are based on two theorems and only two, na mely the existence theorem of Cauchy 



and the theorem of perturbations of Poincare |43]. This is explained in detail in the 



lecture notes of a Chamonix school 0, an updated version of which is in prepara- 



tion 

For discrete equations, this is also the case for all methods but one, the singu- 
larity confinement method, which really seems outside the scope of the theorem of 
Poincare. 

Consider an arbitrary discrete equation (0), also depending on some parame- 
ters a, and let (x, /i, m, a) (X, if, f/, A, e) be an arbitrary perturbation admissible 
by the theorem of Poincare (which excludes any nonanalyticity, like e^^^, for the 
perturbed variables). A necessary condition is that the limit £ = possesses the PP 
(discrete or continuous, this does not matter). 

To illustrate the different methods, let us discretize the equation (PI) by a second 
order scheme, using the rules previously stated 

E = -(II-2m + m)/i-2 + 3Ai(m + m)m + 6A2M^ + 6A3MM + ^ = 0. (69) 
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with = 1, and g an unspecified function of x. 

The test will generate necessary conditions on (A^, g). One must find at least the 

following solution, where the equation has a Lax pair : g = x and A = (2/3, 1/3, 0); 
one may also find the following solutions, isolated by the singularity confinement 

method : ^ = x and A" = (1,0,0) (Ref. [H), T = (1/2,1/4,1/4) (Ref. Q 
eq. (5.5)). 

8.1 Singularity confinement method 

If the field u admits a pole at some point xo in the complex plane x 

u{x) ~ uqx^i X = X — xo 0, Mo 7^ 0, —peM, (70) 

it is generically regular at any point Xq + Xi where Xi is not infinitesimal 

Vxi, |xi| >> : m(xo + Xi) 7^ oo. (71) 

When u satisfies a discrete equation of order N, the implementation of this "con- 
finement condition" consists in requiring the property ( [Til) for + 1 consecutive 
iterates, which generically ensures the property for the next iterates. The polar 
behaviour is then only sensitive during a finite number of iterations. 



8.2 Method of perturbation of the continuum limit 

Defined by an expansion of m as a Taylor series in the lattice stepsize 

+ 00 

X unchanged, h = e, q = e^, u = e"'u^'^\ a = analytic {A,e), (72) 

n=0 

this perturbation generates an infinite sequence of differential equations E^"^^ = 



E = J2 ^"^^"^ (73) 

n=Q 

E(")(x, . . . , «(")) = E(°)(x, m(°))'m(") + i?(")(x, . . . , m("-i)) = 0, n ^7^) 



whose first one n = is the "continuum limit". The next ones n > 1, which are 
linear inhomogeneous, have the same homogeneous part m*-"-* = independent 
of n, defined by the derivative of the equation of the continuum limit, while their 
inhomogeneous part i?*-"-* ( "right-hand side" ) comes at the same time from the non- 
linearities and the discretization. 

This perturbation of the continuum limit is entirely analogous to the pertur- 
bative method of the continuous case, either in its Fuchsian version IT^ or in its 



nonFuchsian one depending on the nature, Fuchsian or nonFuchsian, of the 
linearized equation E^^^ = at a singulier point of u^^\ 

The simplicity of the method is best seen on the Euler scheme for the Bernoulli 
equation ||T^ 

E={u{x + h)-u{x))/h + u{xf = 0, (75) 

i.e. the logistic map of Verhulst, a paradigm of chaotic behaviour which should 
therefore fail the test. Let us expand the terms of (^) according to the perturbation 
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( |72D up to an order in e sufficient to build the ffist equation E^^^ = beyond the 
continuum hmit E^^'^ = 

u = «(o) +^,(1)5 + 0(52) (76) 

= M(°)'+2M(°)n«£ + 0(£=^) (77) 
u{x + h) = u{x)+u'{x)h + {l/2)u"{x)h^ + 0{h^) (78) 



h 

The equations of orders n = and n = 1 



(79) 



= ^(0)' + ^(0)^ = (80) 
^(1) = ^(o)'^(i) + (1/2)^(0)" = 0, = 9, + 2mW. (81) 

have the general solution 

■u*-"'' = X \x = ^~2;o, Xo arbitrary (82) 
m'-^-' = M^lx^^ — Log V^, '?/' = x — Xo, arbitrary, (83) 

and the movable logarithm proves the instability as soon as order n = 1, at the 
Fuchs index i = —1. 

Remark. The only restriction on u'^^'> is not to be what is called a singular solution 
(not obtainable from the general solution by assigning values to the arbitrary data), 
i.e. it can be either the general solution (as above) or a particular one, it can also 
be either global (as above) or local (Laurent series). 

The processing of the example (p9| ) [|T^ isolates three values of A , with g = x. 

The ffist value A = (2/3,1/3,0) (case a = 1 in Ref. [^) corresponds to the 
d-(Pl) (|1|) with a Lax pair found in quantum gravity, the condition is then sufficient. 
The second value (1,0,0) (case a = in Ref. ||2^) corresponds to a candidate d- 
(Pl) with a second order Lax pair [|l^. The third value (1/2, 1/4, 1/4) defines an 
equation equivalent to that for (1,070) under a discrete birational transformation 

mm- , ^ 

Remark. Following Painleve one should in fact search for the "complete 
equation" , i.e. for all the admissible nondominant terms which can be added to the 
candidate d-(Pl) (^) without destroying the PP. According to the already known 
difference or g— difference d-(Pl) candidates (fifteen to our knowledge), the only 
admissible complementary terms seem to be 

/i^Xdiscr(M), h'^udiscT{u'^), {u - u) / {2h) , 

h^X\ h^XdiscT{u^), /iVdiscr(M2), (84) 

in which X is x for difference equations or the suitable exponential function of x for 
g— difference equations. 



8.3 Comparison of the two main methods 

The two main methods which define the discrete Painleve test, namely the sin- 
gularity confinement method and the perturbation of the continuum, happen to find 
the same necessary conditions when applied to a sample of equations : the candidate 
d-(Pl) (|69|) , a candidate discrete Chazy equation of class III not yet fully integrated 
Q, discrete nonlinear Schrodinger equations or various discrete Korteweg-de Vries 



equations Other examples are currently under investigation to detect situations 



where the two methods would produce complementary, not identical results. 
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Since it only relies on the existence of a continuum limit, our method can be 
extended without difficulty to equations in an arbitrary number of dependent or 
independent variables, whether the equations be discrete or mixed continuous and 
discrete. 

The singularity confinement method has also been extended to such situations 



5T|; however, in the case of m discrete independent variables, one must check in 



addition that the result of the iteration is independent of the path followed on the 
m— th dimensional lattice. 

In the case of equations of second degree and higher, our method is unchanged 
since, again, it relies on the continuum limit, for which this technical question is 
settled. In such a case, the confinement method must make at each step a coherent 
choice of determination then compute the confinement condition. 



9 Towards the discrete Painleve and Gambier equa- 
tions 

The task of finding discrete analogues of the fifty canonical equations of Gambier 
d-(Gn) is, at present time, far from being achieved. Let us give here a brief summary 
of the situation and a few lines of conduct to improve it. 

In the continuous case, the PP has been proved either by explicitly linearizing, 
or by integrating with elliptic functions, or by proving the irreducibility and the 
absence of movable critical singularities. This is equivalent to either linearize or find 
a Lax pair. 

In the discrete case, after having performed the discrete Painleve test in order to 
isolate candidates d-(Gn), one must do the same : either linearize or find a discrete 
Lax pair. 

The only sure informations at our disposal are : the fifty continuous (Gn) equa- 
tions of course, the exact discrete elliptic equations. 

To be precise, according to the conjecture of Section ^, one should look for 
d-(Gn) equations satisfying the following criteria. 

1. Each d-(Gn) passes the discrete Painleve test. 

2. Each d-(Gn) has, like (Gn), order two and degree one. 

3. Each d-(Gn) must either be explicitly linearizable, or possess a ffist integral 
identical to an elliptic equation, or possess a matricial Lax pair [A, B, z) ad- 
mitting as continuum limit a Lax pair (L, M, t) of (Gn). This implies an order 
two for these discrete Lax pairs. 

4. The d-(Pn) equations satisfy a confluence cascade admitting as continuum 



limit the cascade of Painleve and Gambier ^ down to the Weierstrass 
level included, see formulae below. 

As an example of some difficulties, consider the (G27) equation in the particular 
case of Ermakov and Pinney 

-uu" + (1/2)m'2 + f{x)u^ - cV2 = 0, c 7^ 0, c constant, (85) 
which is linearizable either by derivation 

-u" + 2fu' + fu = 0, (86) 
or by the singular part transformation 

= c(Log(V'2M))', < - UmiJk = 0, A; = 1, 2. (87) 
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Three discrete candidates for (EB) have been proposed, the first one of degree 



one ^ 

,-2 I ^^ /n\f^ „,\U-'i 



-u{u -2u + u)h-^ + {l/2){u- u){u - u)h-^ + /(u + + m)/4 - c72 4^ 
the second one also of degree one 



-u(u- 2u + u)h ^ + (1/2) (m- u){u - ujU'"^ + fuu- c^/2 

+ {h/2)c{u-2u + u)h-'^ = (89) 



the third one of degree two ||^, |T^ 



(H - M - (2 + h^f/2fuf - (2 + h^f/2f{4:uu + h'^c^)) /{8h^) = 0. (90) 

The third one has been hnearized by discrete analogues of both transformations 
and (^) |T0[ but it has degree two. The first one has been linearized by a 
discrete analogue of (|^) only, but it has two nice features : it obeys the rules and 
has no complementary terms. The second one, to our knowledge, has not yet been 
integrated. 

Before proceeding, let us recall for later use the definition of (P4) 

u" = — + -u^ + Axu"^ + 2(x^ - a)u + - (91) 
2u 2 u 

and the restriction to the subset ((P4), (P2), (PI), p) of the confluence cascade 
from (Pn)(x, M,a;,/?,7,(5) to (Pm)(X, ?7, A, 5, (7, D), m < n, 

4^2: {x,u, a, p) = {eX - e'^/A, e'^/A + e'^U, -£~^/32 - 2A, -£^^^51^^,2) 
2 1 



{x, u, a) = {-66-^^ + e^X, e'^ + eU, 4e-^^), (93) 
(x, u) = {-6-^92/2 + eX, e-'U), (94) 



with e — > 0. 

Let us from now on restrict to discretizations without complementary terms, 
such as, for d-(Pl), equation ( |69D but not equation (|6^. 
For d-(Pl), three candidates ( |69D pass the test : 

T = (2/3, 1/3, 0), (1, 0, 0), (1/2, 1/4, 1/4). Among these, the first one is the only one 
to admit a confluence to the discrete equation of Weierstrass and, since it satisfies 
the other above criteria, it definitely can be called "the (unique) d-(Pl) without 
complementary terms" . 

For d-(P2), this is (H) the good equation, because of both its Lax pair and the 
confluence from d-(P2)(u, x, /i, a) to d-(Pl)(f/, X, if) this confluence is better 



written as 

{x,h,u,a) = {X\X-Qe-^^),X^H,X{e-'^ + U),4^^e-^^), X'"" = e'^ + Hh-^\ (95) 



which proves that its continuum limit is the continuous confluence 

For the four other equations (P3), (P4), (P5), (P6), there is not yet a fully 
satisfactory result, i.e. for each of them at least one of the enumerated points is not 
satisfied. 



Two discrete (P3) candidates have been proposed. The g— (P3) candidate 



has a matricial Lax pair |4T[] of order four, not two, without clear continuum limit. 
The d-(P3) candidate, systeni (25ab) of two equations in two fields in Ref. [^, has 
a fourth order Lax pair [^, |3^, but it also has a subtle drawback. Indeed, after 
elimination of the second field, the discrete equation has degree two, although its 
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continuum limit (P3) has degree one. The second order Lax pair ([58| ) is a starting 
point to remedy this situation. 
The d-(P4) candidate |8 



E = -u{u -2u + u)h-^ + (1 /2) (u -u){u- u)h'^ 

+u^{uu + uu + uu)/2 + {xu + /2){u + u){u + u) - 2au^ + P = 



(96) 



admits a confluence to the d-(P2) and one can even check that this confluence 
is independent of the stepsize and given by (P^. Although this is evidently the good 
d-(P4), its Lax pair is still unknown. 

No d-(P5) candidate is known. The g— (P5) candidate |^ has the correct degree 

_ _ _ 



(one), it admits confluences 
its Lax pair is unknown. 



to both the g— (P3) candidate and the d-(P2), but 



Two discrete (P6) candidates have also been proposed. Both the g— (P6) candi- 
date [0 and the d-(P6) candidate are defined as a system of two equations in 
two fields and, just like for the d-(P3j candidate, the discrete equation in a single 
field has second order but also second degree (as remarked in Ref. ||3^), although 
its continuum limit is (P6) itself. The d-(P6) candidate has not yet a Lax pair, 
but the g— (P6) candidate has one by construction, since Jimbo and Sakai started 
in fact from a g— difference isomonodromic deformation problem in order to obtain 
their system. 



10 Conclusion 

These precise definitions and guidelines may render more systematic the search 
for discrete analogues of differential equations integrable in the sense of Painleve. 

The most fundamental anchor point consists of the discrete elliptic equations 
because they are exact. 

As to the new perturbative method for the discrete Painleve test, its full applica- 
bility (to any number of independent variables, whether discrete or mixed discrete- 
continuous) should make it a quite efficient tool to investigate which discretizations 
of partial differential equations preserve the Painleve property. 
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